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ABSTRACT 



Context. We develop a two-flow model of accretion onto a black hole which incorporates the effect of the local magneto-rotational instability. 
The flow consists of an accretion disk and an accreting corona, and the local dynamo affects the disk/corona mass exchange. 
Aims. The model is aimed to explain the power spectrum density of the sources in their soft, disk-dominated states. 

Methods. The local perturbations of the magnetic field in the disk are described as in King et al. (2004) and Mayer and Pringle (2005), but the 
time-dependent local magnetic field is assumed to aff^ect the local supply of the material to the corona. Since the viscous timescale in the corona 
is much shorter than in the disk, the local perturbations are not smeared in the corona. Simple analytical estimates and full time-dependent 
computations of the disk-corona system are performed. 

Results. The accreting corona model can reproduce the broad power spectra of Soft State X-ray binaries and AGN. The model, however, 
predicts that (i) sources undergoing radiation pressure instability (GRS 1915+ 105) should have systematically steeper power spectra than other 
sources, (ii) AGN should have systematically steeper power spectra than GBH, even if their disks are described using viscosity proportional to 
the gas pressure. More precise measurements of power spectra of Soft State sources are clearly needed. 



Key words, accretion disks ~ X-rays:binaries - galaxies: active 



1. Introduction 

The X-ray emission of the Galactic Black Holes (GBH) and 
active galactic nuclei (AGN) is strongly variable. Aperiodic 
character of this variability (GBH: Terrell 119721 for a review 
see van der Klis 1989 Remillard & McClintock |2006l AGN: 
McHardy & Czerny [19871 Lawrence et al. [T9871 for a re- 
view, see Mushotzky et al. '1993) makes the interpretation of 
the physical nature of the phenomenon rather difficult. On the 
other hand, several important properties of this variability were 
found: (i) the shape of the power spectrum depends signifi- 
cantly on the luminosity state of the system (e.g. Pottschmidt 
et al. 120031 McHardy et al. l2004l i: (ii) in Soft State sources the 
slope of the power spectrum is close to 1 in a broad frequency 
band, only steepening off at the highest frequencies which 
means that a broad range of frequencies contribute roughly 
equal power to the overall variability (e.g. Revnivtsev et al. 
I2OOOI 1: (iii) in Hard State sources the power spectrum is more 
complex (although the energy spectrum is much simpler, be- 
ing roughly of a power law shape with energy index ~ 0.9), 
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well modeled by a few Lorentzians with a power law compo- 
nent also present (e.g. Pottschmidt et al. 120031 1: (iv) the sim- 
plest shot noise cannot account for the brightest flares unless 
the flares come in bunches or avalanches (Canizares & Oda 
[T9771 Negoro et al. [T995l Poutanen & Fabian [T999l l; (v) the un- 
normalized excess variance scales with the flux (see e.g. Uttley 
& McHardv l2001!) . The Hard State sources are possibly filtered 
versions of the Soft State variability, therefore the understand- 
ing of the Soft State variability is of key importance. 

Several models of X-ray variability were considered: rotat- 
ing hot spots at the accretion disk surface (Bath et al. 1 19741 
Abramowicz et al. I1992I I. magnetic flares above the accretion 
disk (Galeev et al. 1 19791 Poutanen & Fabian 1 1999! Goosmann 
et al. I2OO6I 1. variable obscuration (e.g. Boiler et al. 119971 
McKernan & Yaqoob [19981 Abrassart & Czerny |2000l Karas 
et al. lIOOO l. local turbulence (Nowak & Wagoner [19951 ). But 
the key problem in these models remained: why long timescale 
variability, plausibly related to the processes in the outer parts 
of the disk contributes so significantly to the overall variability 
although the energy dissipated in the outer parts of the disk is 
very low? 
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A promising solution was suggested by Lyubarskii (|1997| l 
who postulated that the basic variability mechanism operates 
locally at each disk radius, as fluctuations of the viscosity pa- 
rameter a caused by the dynamo mechanism, and these per- 
turbations propagate inwards, finally modulating the accretion 
rate in the innermost part of the disk. Therefore, these modu- 
lations, containing the memory of the whole disk perturbations 
are finally seen as X-ray emission variability. Such a mecha- 
nism in a natural way explains the slope of the power spectrum 
in a disk extending down to the marginally stable orbit. 

However, the propagation of the perturbations, as envi- 
sioned in Lyubarskii model, is not a straightforward issue. In 
the original model, local dynamo effect was expected to gen- 
erate local fluctuations in the accretion rate, and these fluc- 
tuations were expected to propagate inwards. The problem is 
that these fluctuations propagate in a viscous time, and they 
are smeared unless the local dynamo timescale is longer than 
the local viscous timescale, as pointed out by Churazov et al. 
( l200n see also Kotov et al. 2001 Gleissner et al. l2004] l. These 
problems were overcome by combining the local magnetic dy- 
namo with magnetic field difiiision and wind/jest outflow (King 
et al. 120041 Mayer et al. l2006] l. 

In the present paper we propose to accommodate the ba- 
sic picture of Lyubarskii within the frame of the accreting 
corona model. In our model the variations due to the local dy- 
namo modify the accretion rate within the corona. The disk 
perturbations propagate in its viscous timescale and therefore 
are smeared, but the viscous timescale in the corona is much 
shorter than in the disk. This allows the coronal perturbations 
to be observed. The model accounts for the correct power spec- 
trum of the X-ray lightcurves of sources in the soft states as 
well as for the fact that in GBH disk flux seems to remain al- 
most constant and all the variability is likely coming from the 
corona (Churazov et al. 1200 11 1. 

The structure of the paper is the following: the simple semi- 
analytical model is discussed in Sect. |2] the results based on 
numerical model of the accreting disk - accreting corona sys- 
tem are described in Sect.|3] the discussion in the context of the 
observational data is presented in Sect. IH and we conclude in 
Sect.|5] 



2. 1 . Power law trends 

We will start with simple analytical estimates of the power 
spectrum of the X-ray lightcurve. This is possible if the de- 
pendence of the process on the disk radius is in the form of a 
power law, and the effect of the inner and outer radius can be 
neglected. There are three key quantities built into the model. 

First, we assume that the variability is caused by some pro- 
cess which causes the coherent changes in the accretion disk 
properties which span certain range of radii, Ar. The size of the 
region is likely to depend on the location of the perturbation. 
We assume 

ArocA. (1) 

Next, the characteristic timescale of the process depends on the 
location of the perturbation as 

tocr^'. (2) 

Finally, the amplitude of the ejection rate from the unit disk 
surface element, in-, is also likely to depend on the location. 
We assume 

m, oc r"' . (3) 

In such a model, the ejection of mass in a single event has 
an amplitude A oc ih,rAr, if the event happens in a whole ring, 
and the slope of the power spectrum describing the variability 
is then given by 

P(/)./-.; 1,111^. (4) 

If we want to reproduce the slope p - 1 characteristic of the 
Soft State X-ray sources, we should choose models satisfying 
the condition 

2a r + = -3, (5) 

and the condition is independent from y, , i.e. the way how 
the timescales scale down with the radius is unimportant, as 
stressed by Lyubarskii (1997). 

This analysis neglects the propagation time in the corona, 
and the effect of the boundary conditions. 



2. Semi-analytical model of accreting corona 

We consider here the case of an accretion disk extending down 
to the marginally stable orbit, surrounded by the hot optically 
thin corona. We postulate that the corona itself can also trans- 
port the angular momentum and can exchange mass and an- 
gular momentum with the underlying optically thick disk. The 
model was proposed in the context of disk evaporation in cat- 
aclysmic variables (Meyer & Mever-Hoffmeister 119941 and 
subsequently developed in a number of papers (Zycki et al. 
[19951 Witt et al. [T9971 Janiuk et al. l2000l Rozahska & Czerny 
120001 Meyer-Hofmesister & Meyer|200T] Misra & Taam|2002l 
Meyer-Hofmeister et al. i2005| . 

Here we simplify the model by assuming no condensation 
of the coronal material onto the disk, i.e. the disk material is 
continuously ejected into the corona, according to the adopted 
laws. 



2.2. Physical semi-analytical model 

We now develop a more complex, semi-analytical approach 
which nevertheless captures the essential physics of the accret- 
ing corona model. We assume that the dynamo mechanism is 
operating within the disk and we model the stochastic aspect of 
this process as suggested by King et al. (I2004I I. 

We construct the discrete radial grid corresponding to 
the radial size of the perturbation. We consider two cases: ei- 
ther Ar - H (case a), where h is the local disk thickness taken 
from a stationary model, or Ar = r (case b), which represents 
the scaling with the thickness of the accreting corona. Here 
Ar corresponds to the difference between the consecutive radii, 
rk+\ - fk- 

At each disk radius r^ we create a time series assuming an 
underlying Markoff process 

M„+i = -a\Un + en, (6) 
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where e„ is a random variable between - VS and V3, which has 
a uniform distribution, zero mean and the dispersion equal to 
unity. We assume a\ - 0.5 after King et al. (120041) and Mayer 
& Pringle (l2006l l. 

The time-step unit at each radius is proportional to the dy- 
namical (orbital) time at a given disk radius: 



Table 1. Semi-analytical models 



* mag 



GM' 



(7) 



where M is the black hole mass, and the factor is the pa- 
rameter of the model. The parameter kd can be considerably 
larger than 1 since the magnetic processes are slowed down by 
macroscopic resistivity (see e.g. Tout & Pringle [T992[ Stone et 
al. .l996. Komissarov et al. ,2006b . 

The Markoff chain formed at each radius describes the time 
evolution of the magnetic field 



B - u„B„ 



(8) 



where B„ax(r) is the variability amplitude of the magnetic field. 
This local normalization factor B,nax{r) is determined by the 
effective viscosity parameter a of Shakura & Sunyaev 1 19731 ) 



Bi 



4-n 



= aP, 



(9) 



where P is the total (gas plus radiation) pressure within the disk 
interior. 

Now we estimate the amount of outflow of the disk mate- 
rial into the corona, m,. We postulate that the local disk mag- 
netic field is responsible for the outflow of the material from 
the disk surface to the corona. This outflow is a combined re- 
sult of the magnetic field buoyancy and the numerous possi- 
ble mechanisms of the coronal heating (e.g. magnetic recon- 
nection, Galeev et al. 119791 dissipation of Alfven waves in 
the corona, Kuperus & Ionson [T985; acoustic heating, coupled 
with convection, Bisnovatyi-Kogan & Blinnikov 119771 verti- 
cal Poynting flux, Merloni [20031 electron conduction, Meyer 
& Mever-Hofmeisterl 1 9941 Rozanska & Czerny |200 0', ion irra- 
diation, Deufel & Spruit l2000] l leading to the disk evaporation. 
Since the exact mechanism is under discussion, we adopt three 
viable parameterizations of the evaporation process: 



8;r kT„, 
B- Mm,' 



1 

/('•) 



(case A) 
(case B) 
(case C) 



(10) 



Here we assume that the coronal (ion) temperature is equal to 
the virial temperature, va is the Alfven speed (v^ - P/p, where 
p is the density of the disk), and f(r) is the boundary condition 
factor (equal to 1 - y]?>Rschwlr in the Newtonian case). The 
dependence of the pressure P, and the density p, is taken from 
the analytical formulae of Shakura & Sunyaev (1973), either 
for the gas pressure dominated region (for opacity dominated 
by electron scattering), or for the radiation pressure dominated 
region. 



Model 


m- 


cell 


SS model 


slope 


Aa-rad 


A 


a 


radiation pressure 


0.66 (curv) 


Ba-rad 


B 


a 


radiation pressure 


2.42 


Ca-rad 


C 


a 


radiation pressure 


0.48 


Ab-rad 


A 


b 


radiation pressure 


1.16 (curv) 


Bb-rad 


B 


b 


radiation pressure 


2.99 


Cb-rad 


C 


b 


radiation pressure 


1.00 


Aa-gas 


A 


a 


gas pressure 


1.19 (curv) 


Ba-gas 


B 


a 


gas pressure 


1.75 


Ca-gas 


C 


a 


gas pressure 


1.03 


Ab-gas 


A 


b 


gas pressure 


1.15 (curv) 


Bb-gas 


B 


b 


gas pressure 


1.71 


Cb-gas 


C 


b 


gas pressure 


1.00 



Evaporation parameterized as in Eq.[TO](cases A, B, and C); cells a and 
b refer to the size of the magnetic field perturbation, Ar, equal either to 
the disk thickness h(r) or to the current radius, r, respectively. The disk 
structure is described as in Shakura & Sunyaev (1973), either in the 
radiation pressure dominated region or in the gas pressure dominated 
region, but with the electron scattering opacity. The slope of the PSD 
is measured in the low frequency part. The note (curv) means that a 
certain level of curvature is seen in the long wavelength part of the 
PSD. 

We further assume that the coronal inflow proceeds in the 
viscous timescale of the corona, which is taken as 



1 



r3 

GM' 



(11) 



where acor is the viscosity parameter characterizing the corona. 
We allow for the dispersion in a flow of each clump of the coro- 
nal material by a fraction of the inflow timescale, 6. We assume 
that the coronal material does not settle down onto a disk, so 
the effective accretion rate at the inner disk radius consists of 
the coronal clumps ejected from the disk at each radius under 
consideration. 

Finally, the X-ray luminosity of the flow is assumed to be 
proportional to the coronal accretion flow at the inner disk. The 
lightcurve is analyzed using the standard FFT package and the 
resulting power spectrum is binned for clarity. 

2.3. Results 

Since we consider three options for the disk evaporation effi- 
ciency, two options for the size of the elementary cell, and two 
cases of the disk structure, we finally have a family of 12 mod- 
els. We summarize them in Table[T] 

We can understand the results by referring to the power law 
dependencies given by Eq. |4] For case (A) and case (C), the 
important factor is the energy flux provided with the Alfven 
speed to the corona, B^/(87r)v^. Since we assume the scaling 
of the magnetic field with the pressure given by Eq. |9l it can 
be easily shown with the use of the hydrostatic equation for the 
disk structure that for a stationary Shakura-Sunyaev disk the 
amplitude of this factor reduces to 



—VA oc F{r) oc — — , 



(12) 
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Fig. 1. Power spectrum density for the radiation-pressure dom- 
inated and gas-pressure dominated disk if the dynamo cells 
scale with the disk thickness (models Aa-rad and Aa-gas from 
Table[B. 



Fig. 2. Power spectrum density for the radiation-pressure dom- 
inated (dotted histogram) and gas-pressure dominated (solid 
histogram) disk if the dynamo cells scale with the disk radius. 
Random location of cells (models Ab-rad and Ab-gas from 
TablelB- 



where F(r) is the radiation flux of a stationary disk, and f(r) 
is the boundary condition factor. Therefore, this value does not 
depend on the disk model, i.e. whether the disk is dominated 
by the radiation pressure or the gas pressure. The radial shape 
of the dependence is then given by oc r^^ fir) for the 
case (A) and it is just m, oc r^^ for the case (C). Therefore, 
case (C) is characterized roughly by a power law dependence 
with Ur - -2, and case (A) shows a departure from a power 
law behaviour due to the f{r) factor. This means that the inner 
boundary condition directly affects the shape of the PSD. 

The resulting power spectrum strongly depends on the disk 
thickness, if the dynamo cells scale down with it. As an illus- 
tration, we show in Fig. [1] the two power spectra: one for ra- 
diation pressure dominated disk and one for the gas pressure 
dominated disk. This means that the power spectra of AGN 
with disks likely to be dominated by radiation pressure should 
show systematically smaller values of the slope of the power 
spectrum, p. 

However, if we assume that the magnetic cell is large and 
equal to the current disk radius (Ar - r), the power spectrum 
is the same for the gas pressure and for the radiation pressure 
disk model (see Fig. |2|i. There are only a few cells, and if the 
cells were always located at the fixed disk radii, the power spec- 
trum would show clear regular peaks. We may overcome this 
numerical problem by means of randomization of the flare lo- 
cation around the predicted location or by the assumption of 
the thinner cells (e.g. dr - 0.05r). This smears this effect out 
and the resulting power spectrum, shown in the Figure, does not 
show specific features. The resulting power spectrum is slightly 
steeper than p - I. 



Modification of the energy and mass transfer between the 
disk and the corona also affects the PSD slope significantly. 
Models based on Poynting flux (B-type) are always very steep, 
dominated by the longest timescales. C-models, on the other 
hand, show no sign of curvature at the long wavelengths, and 
they are always flatter. Three of these models have the PSD 
slope p very close to 1. 

3. Time evolution of the disk with accreting corona 

3.1. Model 

The model of an accretion disk exchanging the mass with the 
hot corona is developed from the model described and studied 
in detail in Janiuk & Czerny (2005). The basic equations of the 
(vertically averaged) accretion disk model are: (i) the viscous 
heating 



(ii) the radiative cooling: 



rad 



(13) 



(14) 



where we adopt the electron scattering opacity k - 0.34 cm^/g; 

(iii) the advective cooling (e.g. Paczyhski & Bisnovatyi- 
Kogan il981t Muchotrzeb & Paczyhski |T982j Abramowicz et 
al. [T988]) 

2rf ^adv „ 



adv 



3pGM 



(15) 
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and ^adv ~ 1 -O. 

Similarly to the semi-analytical calculations presented in 
Section|2l we differentiate between two cases: gas pressure and 
radiation pressure dominated disks. Therefore, we adopt two 
kinds of heating prescriptions in Eq. [T3] P - Pgas and P = 
f tot- In the latter case the total pressure consists of the gas and 
radiation pressure, P - -j^pT + ^aT'^, and we use a modified 
viscosity parameterization 



(l+^/^o) 
(1 + (^/^o)2 



(16) 



where ^ = PmdlP^as and - 0.01 and = 8.0 (Nayakshin 
et al. l2000l l. This is because the constant viscosity parameteri- 
zation, appropriate for the gas pressure dominated disk, would 
lead to much too large amplitudes of the disk outbursts due to 
the thermal-viscous instability induced by the radiation pres- 
sure. 

The angular velocity of the disc is assumed to be Keplerian, 
Q. - -ijGM/r^, and the sound speed is Cj = -\JP/p = ^H. 
A non-rotating, Schwarzschild black hole is assumed and the 
inner radius of the disc is always at 3 /?schw The outer radius is 
equal to 300 ^schw, and at this radius we give a constant mass 
inflow, parameterized by the external accretion rate Mext. The 
mass of the black hole is assumed to be IOMq or 10^ Mq. 

To start the disk evolution, we calculate the initial quasi- 
steady configuration from the energy balance: Ftot - Q^i^c ~ 
2adv 2iad' where the total energy flux Ftot dissipated within 
the disc at a radius r is calculated as: 



3GMM 
F,„, = — — —f(r) 



and 



/('•) = 1 



(17) 



(18) 



is the boundary condition in the Newtonian potential. 

The corona above the disk is thick, //coi - r, and the 
pressure in the corona is its gas pressure due to ions: Pcoi = 
^Pcoi-T^cor- The coronal temperature is equal to the virial tem- 
perature, T.or = ^ ^ . 

The initial configuration of the corona is computed under 
the assumption that the corona has a uniform surface density 
and its optical depth is constant: 



0.1 



(19) 



The treatment of the mass exchange in the current model is 
different from the one adopted in the previous paper (Janiuk & 
Czerny 2005), since now we incorporate the stochastic charac- 
ter of the dynamo responsible for the coronal heating. The mass 
exchange is proportional to the locally generated magnetic flux, 
divided by the energy change per particle. Therefore the rate of 
the local mass exchange in vertical direction, in units of g s"' 
cm"-, is given by Eq.[TOl where f(r) is the boundary condition 
given by Eq. [18] va - ■\/aQ.H is the Alfven velocity and the 
poloidal magnetic field component is calculated from 



For generality, we introduced here the parameter f3s < 1 .0. It is 
adopted here to account for the possible diff'erence between the 
average magnetic field in the disk and the larger scale magnetic 
field likely to be more important for the disk evaporation (e.g. 
Lovelace, Romanova & Newman 119941 King et al. I2004I I. In 
the previous section we assumed always /3s = 1, for simplicity. 

The random changes of the magnetic flux are computed as 
follows. First, we assume certain law for the radial distribu- 
tion of the magnetic cells (we consider Ar - H{r) - case a, 
Ar = yfr - case b, and Ar = r - case c). Therefore we impose 
a separate radial grid to follow the stochastic changes of the 
B field in these boxes. Next, we assume that the timescale of 
these stochastic fluctuations is proportional to the dynamical 
timescale: 



Td - Kd/Q., 



(21) 



with Kd = 100. Since our disc is geometrically thin, with a ratio 
H/r as small as ^ 0.01 in some of the small mass models, and 
even less in the supermassive black hole models, the number of 
points in the magnetic grid varied from ~ 150 to ~ 2250. 

Having computed the initial disc and corona state we solve 
the equation of mass and angular momentum conservation: 



^ = --(3ri/2 (^i/2yS))-mz 
ot r or or 

and the energy equation: 

dT dT T 4-3B dl. dl, 

+ v.- ( + Vv — ) 

dt dr E12-10.5yS 5f dr 



(22) 



(23) 



1 



+ - 



PH 12- 10.5/3 

Here S = Hp is the surface density in the disc, - 
j^^(vSr'^^) is its radial velocity, and v = (2Pa)/(3pQ) is 
the kinematic viscosity. The heating term is given by Equation 
[T3] and the cooling term Q is given by Equation [14] while 
the advection is included in the energy equation via the radial 
derivatives. 

The evolution of the coronal density is also given by mass 
and angular momentum conservation: 



^ = i|-(3r'/2|.(,i/2, 
ot r or or 



tScoi)) + OTz 



The radial velocity in the corona is vj:'^ 



(24) 



fAI2 ^^-V ^coi-'-'cor' } 



(20) 



and viscosity is Vcm- = (2Pcor<3'coi-)/(3pcorf^) where acor = 0.01. 
The coronal temperature is always equal to the virial tempera- 
ture and does not vary with time. 

We solve the above set of three time-dependent equations 
using the variables y - 2r'^^ and E = at the fixed radial 
grid, equally spaced in y, and the number of radial zones is 
set to A^i - 100. For the time evolution we use the fourth or- 
der Runge-Kutta method and the Adams-Moulton predictor- 
corrector method. We choose the no-torque inner boundary 
condition, Ein = Tin = for the disc. The outer boundary of 
the disc is parameterized by an external accretion rate M^w 
The boundary condition in the corona is given by Equation[T9] 

with rcor(/?out) = Tcor(^m) =0.1. 

The local mass exchange rate depends on time via the 
stochastic changes of the magnetic flux. These are calculated 
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Table 2. Numerical models 



Model Mass [Mq] Instab. m- cell size Heating slope 



1 


10 


+ A 


c 




1.25 


2 


10 


+ B 


c 


Aol 


1.90 


3 


10 


+ c 


c 


^tol 


0.99 


4 


10 


+ A 


a 




1.65 


5 


10 


+ B 


a 


^tot 


2.75 


6 


10 


+ C 


a 


^tot 


0.89 


7 


10 


A 


c 


Aol 


1.13 


8 


10 


B 


c 


Aot 


1.77 


9 


10 


C 


c 


^tot 


0.85 


10 


10 


A 


a 


^tot 


1.60 


11 


10 


B 


a 


^tot 


2.25 


12 


10 


C 


a 


^tot 


1.37 


13 


10 


A 


c 




1.02 


14 


10 


B 


c 


P 

' gas 


1.71 


15 


10 


C 


c 


^'ga. 


0.38 


16 


10 


A 


a 


^'gas 


0.94 


17 


10 


B 


a 


^'gas 


1.88 


18 


10 


C 


a 


^gas 


0.58 


19 


10** 


A 


c 


P 

' gas 


1.02 


20 


10** 


B 


c 


^gas 


2.96 


21 


10" 


C 


c 


P 

' gas 


0.82 


22 


lO'' 


A 


a 


^'gas 


0.95 


23 


10** 


B 


a 


/"gas 


1.92 


24 


10" 


C 


a 


P 

' gas 


0.78 


25 


10^ 


C 


b 


^gas 


0.93 


26 


10 


C 


b 


^gas 


0.76 


27 


10» 


A 


b 


^'gas 


1.59 


28 


10 


A 


b 


^'gas 


0.86 


29 


10** 


B 


b 


^gas 


2.97 


30 


10 


B 


b 


/'gas 


1.72 



The prescriptions for disk evaporation: A, B, C, as given by Eg. 1101 
The cell size (a) is Ar oc H, the size (b) is Ar, and the cell size (c) 
oc s/r,. The viscous heating is proportional either to the total pressure 
(with modified viscosity law) or to the gas pressure. Whenever the disk 
is unstable due to radiation pressure, and the cell size is proportional 
to the disk thickness, the H(r) profile is taken from the hot state of the 
disk. The slope of the PSD is measured in the low frequency part. 



on the magnetic cell's grid with s> A^, . Whenever the time- 
step exceeds the local magnetic timescale (Eq. l2Tl i. the new 
value of a random number is created, using the MarkofF series 
according to Eq. |6] Because of a different number of points in 
the two grids, A^ni and A^i , the local mass exchange rate on the 
basic grid A^i- can change 0, 1 or a few times in one time-step. 
Therefore we simplify this problem, and if the magnetic cells 
are distributed much denser than the basic radial grid points, 
and ih^ would change a few times, we take the average of the 
random numbers generated for one point in the same time- 
step. In particular, this is the case in the innermost parts of the 
accretion disk. 

3.2. Results 

In Table |2] we list the calculated numerical models and their 
resulting power spectral slopes. 




log f [Hz] 



Fig. 3. Thick histogram: PSD from the full evolution of 
a disc-corona system with the mass exchange (model 12 
from Table |2]l ; dotted histograms: power spectra from semi- 
analytical scheme (lower histogram: model Ca-gas, upper his- 
togram: model Ca-rad); thin lower histogram: power spectrum 
from semi-analytical model but with numerical prescription for 
the disk thickness and the total pressure distribution from a 
stationary numerical model; dotted line is representation by a 
power of the slope 1 .2 and the high energy cut-off 10 Hz. 



The slopes were calculated by means of fitting a power law 
function to the numerically computed PSD, in the frequency 
range log / = (-1.8; -0.8) Hz and (-8.0; -7.0) Hz for Galactic 
black holes and AGN, respectively. The calculations were per- 
formed for stellar mass accreting black hole with M - IOMq 
(gas and radiation pressure dominated disks) and supermassive 
black hole with M - 10*^ Mq (only gas pressure dominated 
disk). We adopt the parameterizations of the evaporation pro- 
cess as listed in Eq. [10] (cases A, B and C), with different sizes 
of the magnetic cell: (a) Ar oc H, or (b) Ar oc r. We also cal- 
culate the models with size of magnetic cell proportional to the 
square root of radius: (c) Ar oc ^r, which results from our ra- 
dial grid. 

The power spectrum resulting from the fuU time evolution 
study of the disk/corona system systematically differs from the 
simplified semi-analytical model. We can see from Fig. [3] that 
the analytical model with the same physical assumptions for the 
disk evaporation rate gives much flatter power spectrum (dotted 
histograms) than the numerical model. 

One of the reasons is that the disk properties in the numer- 
ical model differ from the analytical prescription. Numerically 
calculated disk under assumption of the total pressure is not 
fully radiation-pressure dominated, therefore neither the total 
pressure, nor the disk thickness in analytical model actually fol- 
low the numerical results. As an example, we plot in Fig.|4]the 
ratio of the disk thickness in the numerical model and in the an- 
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Fig. 4. The ratio of the disk thickness to the radius for numer- 
ical model (continuous line) and analytical radiation pressure 
and gas pressure models (dashed lines) shown in Fig. [3] 



Fig. 5. Power spectra calculated from numerical models with 
the prescription A for mass evaporation rate (cf. Eq. [TO] and 
Tab. 121 the labels on separate curves refer to the number of the 
model as given in the Table). 



alytical gas-dominated as well as radiation-pressure dominated 
disk model. Since in the size of the cells in (a) type models for 
a galactic object is rising faster than linearly with the radius, 
the variability at the outer edge gains systematically additional 
power with respect to the variations in the inner part of the 
disk and the power spectrum resulting from numerical model 
is steeper than from the analytical one. Therefore we also cal- 
culated a semi-analytical model but with the disk thickness and 
the total pressure distribution taken from a numerical station- 
ary model. The result is shown in Fig. [3] as a thin histogram. 
The PSD slope of this model is p = 1.13, i.e. in the numerical 
model the PSD is slightly steeper due to the exact computation 
of the propagation of perturbations through the corona. 

The numerical computations include one more extremely 
important effect: they take into account the non-stationarity 
of some disk models. It is well known that the heating pro- 
portional to the total pressure leads to the disk instabilities 
and luminosity outbursts, if only the external accretion rate is 
large enough to make the radiation pressure dominate (Pringle, 
Rees & Pacholczvk |1973l Shakura & Sunyaev [19761 see also 
Szuszkiewicz & Miller [19971 [T9981 Nayakshin et. al. |2000l 
Janiuk, Czerny & Siemiginowska l2000l l20D2i Mayer & Pringle 
I2OO6I 1. In our models the instability switches on for rii > 0.03, 
for M - IOMq black hole mass. In this case the disk thick- 
ness strongly fluctuates in the unstable region. It means that 
the cell size also strongly fluctuates, however numerically it is 
too much time-consuming to follow this evolution at each time 
step. Therefore, when we adopt the cell size as proportional to 
the disk thickness, the H(r) profile is taken from the hot state 
solution. 

As in the semi-analytical models, the predicted variability 
significantly depends on the adopted evaporation law. In Figure 



|5lwe plot the power spectra calculated from numerical models 
with the prescription A (see Eq. [TOl l. The total pressure mod- 
els "A" with the cell size proportional to (01, 07) have the 
power spectra substantially flatter than the models with cell's 
of sizes proportional to the disk height (04, 10), respectively 
for other parameters kept the same, and the unstable models 
(01, 04) are steeper than their stable equivalents (07,10). Gas 
pressure models (13 and 16) have flatter power spectra than the 
total pressure models. 

In Figure |6l we plot the power spectra calculated from nu- 
merical models with the prescription B for mass evaporation 
rate. All power spectra are very steep, as in analytical solutions, 
independently from other assumptions. 

In Figure |7l we plot the power spectra calculated from nu- 
merical models with the prescription C for mass evaporation 
rate. This prescription gives the flattest power spectra. The 
trend of spectrum flattening for the cell's size proportional to 
■\/r is present in stable models. The dependence on the heating 
term (Ptot vs. f gas) is similar to the case of models "A". 

We also calculated the disk-corona system evolution and 
power spectra for the case of accreting supermassive black hole 
in AGN. We took the black hole mass of lO^M© and an accre- 
tion rate of m - 0.35 in Eddington units. The heating of the 
gas was always assumed to be proportional to the gas pres- 
sure, since the radiation pressure instability in the AGN disks 
results in much stronger changes in disk density and thickness 
between the hot and cold states. Therefore, the evolution of an 
unstable disk is much more complex from the computational 
point of view, i.e. it requires a substantial fine-tuning of the 
time step, and at the moment it is beyond the scope of our com- 
puter facilities. 
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Fig. 6. The same as in Fig|5] but for the prescription B for mass 
evaporation rate. 
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Fig. 7. The same as in Fig|5] but for the prescription C for mass 
evaporation rate. 



The Figure |8] shows the resulting power spectrum density 
for 6 models with various prescriptions for the mass evapora- 
tion rate and magnetic cell size. 

The influence of the parameterization of magnetic cell's 
size (a, b or c) on the PSD slope is less crucial than the role of 
mass evaporation parameterization (A, Z? or C). In general, the 
cell's size Ar proportional to yfr (case c) gives either similar 
or slightly flatter slopes than Ar oc r (case b), while the flat- 
test slopes are produced in the models with cell's sizes Ar cc H 



Fig. 8. Power spectra calculated from numerical models for su- 
permassive black hole in AGN, with the prescriptions A, B, and 
C for mass evaporation rate (cf. Eg. [TOland Tab. |2]i; the labels 
on separate curves refer to the number of the model as given in 
the Table. 



(case a). The slope of the model 27 is exceptionally steep for 
a class A model but this is connected with the frequency lo- 
cation of the slope determination; a small shift towards lower 
frequencies brings this model back to the general systematic 
trends. 

The exemplary lightcurves of the disk and the corona of 
Galactic black hole binaries are shown in Figures |9] [10] and 
nn In Fig. |9]we show the results for the unstable disk evolu- 
tion (04 in Table |2]i, with a large mean accretion rate and disk 
heating proportional to the total pressure. This model results in 
strong disk luminosity outbursts. The coronal luminosity fol- 
lows these long outbursts, and the short variability fluctuations 
due to the magnetic dynamo is superimposed on the long term 
evolution. The magnetic dynamo parameter, was assumed 
to be 0.5, and nevertheless the coronal fluctuations had a sub- 
stantial amplitude. 

In Figs. [TO]and [TT]we show the lightcurves resulting from 
a stable disk model (13 and 18 in Table |2|, with a small mean 
accretion rate and disk heating proportional to the gas pressure. 
The models differ with respect to the prescription for the mass 
evaporation rate and cell size (A-c and C-c, respectively). The 
disk emission is practically constant but the coronal emission 
shows modulation due to the dynamo process. Here the mag- 
netic parameter, /3s, had to be taken large, in order to produce 
significant coronal variability. 

In Figure [12] we show the lightcurves for the AGN evolu- 
tion (model 24 in the Table). All the models calculated for AGN 
disks assumed the gas pressure domination. The key element of 
the models was the prescription for the mass evaporation rate 
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Fig. 9. Exemplary fragment of the disk and corona lightcurves 
obtained from the evolutionary model of an unstable disk (heat- 
ing proportional to Ptou model 04 in Table |2]i. The model pa- 
rameters: M = IOMq, riiout - 0.26, ySs = 0.5. The total duration 
of the simulation is 46 ks. 



Fig. 10. Exemplary fragment of the disk and corona lightcurves 
obtained from the evolutionary model of a stable disk (heating 
proportional to Pgas; model 13 in Table |2]i. The model parame- 
ters: M - IOM0, mout = 0.026, ySs - 1-0. The total duration of 
the simulation is 48 ks. 



and the cell's size: the model "C-a" gave a significant variabil- 
ity in short timescales. 

4. Discussion 

Power spectrum density (PSD) of the X-ray emission of 
Galactic black holes and active galactic nuclei is usually well 
approximated by a single power law with a slope ~ 1, turn- 
ing down at high frequencies, if the source is in the Soft State. 
This means that a broad range of frequencies contribute to the 
lightcurve, although most of the energy is dissipated close to 
the black hole, where the highest frequencies should dominate. 
In the present paper we propose a two-layer accretion flow - an 
accretion disk with an accreting corona - as a solution to the 
problem. 

The accreting corona model was considered before in a 
number of papers (Zycki et al. 1995 Witt et al. 119971 Janiuk 
et al. i2000i Rozariska & Czerny ,2000 , Meyer-Hofmesister & 
Meyer 120011 Misra & Taam 1^021 Meyer-Hofmeister et al. 
I2OO5I 1. The model has several known advantages. For exam- 
ple, it can explain the phenomenon known as the hysteresis ef- 
fect, i.e. for a given source luminosity the spectra are softer 
when the count rate is decreasing, and harder when the count 
rate is increasing (Smith et al. .2002. Meyer-Hofmeister et al. 
120051 Pottschmidt et al. I2006I I. as well as the time delays of 
the hard X-ray band with respect to the soft X-ray band in 
GRS 1915H-105 (Janiuk & Czerny l2005Tl. 

In this paper we considered the propagation of perturba- 
tions in the accretion flow created locally by the magneto- 
rotational instability. We followed the general approach of 
King et al. 120041 and Mayer & Pringle l2006l but we incorpo- 
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Fig. 11. The same as in Fig. [TO]but for the model 18 in Table 
|2](Pgas heating; accretion rate ih - 0.26, /3s - 0.5 ). The total 
duration of the simulation is 45 ks. 

rated their scheme into our two-layer flow. Since the viscous 
timescale in the corona is short, local perturbations affecting 
the coronal accretion rate propagate without significant smear- 
ing. Therefore, the two-flow model offers also in a natural way 
a physical realization of the suggestion made by Lyubarskii 
(119971 ) that the broad PSD is produced by the inward propa- 
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Fig. 12. Exemplary fragment of the disk and corona lightcurves 
obtained from the evolutionary model of a stable AGN disk 
(heating proportional to Pgas; model 24 in Table |2]i. The model 
parameters: M - 10**Mo, mout - 0.35, j3s - 0.1. The total 
duration of the simulation is 415 yrs. 




2x108 4x10" 6x10" BxlO" 

r [cml 



Fig. 13. The thickness of the disk as a function of radius in the 
unstable model of GBH (f tot heating). The solid and dashed 
lines show the H(r) profile in the hot and cold states, respec- 
tively. The black hole mass is IOMq and accretion rate is 
m = 0.26. 



gation of the perturbations created at various disk radii, so the 
inner dissipation region preserves the memory of all fluctua- 
tions. Propagation through the disk, considered originally by 
Lyubarskii (1997), and later explored by King et al. (120041 1 and 
Mayer & Pringle (2006). requires rather long local magnetic 
timescale; otherwise the fluctuations are washed out in course 
of propagation. In the case of disk/corona coupling the mag- 
netic field fluctuations result in modulations of the coronal flow, 
and the viscous timescale in the corona is much shorter than in 
the disk, so the conditions to preserve the fluctuations are easier 
to satisfy. 

The PSD spectra based on this model were calculated 
numerically, but we also provided a semi-analytical insight. 
Below, we discuss the theoretical predictions of our model, and 
then we compare them with the observational data for the Soft 
State X-ray binaries and Narrow Line Seyfert 1 galaxies. 

4.1. Model predictions 

In the modeling, we use a couple of different scaling laws for 
the magnetic dynamo cell: the cell size is either assumed pro- 
portional to the radius of the disc, or its square root, or to the 
disc thickness. This last option is widely believed to be the most 
physical one since the magneto-rotational instability operates 
in a scale of the disk thickness, or a fixed fraction of it, given ef- 
fectively by the parameter a (e.g. Tout & Pringle [T992| Hawley 
et al. [T995l Balbus & Hawley [T998] l. On the other hand, larger 
scale fields can be created either by inverse cascade mechanism 
and reconnection (Livio ,1996. Tout & Pringle . 1992. King et al. 
2004), or by field generation in the corona. 



When we assume that the characteristic dynamo cell is 
determined by the disk height, the resulting power spectrum 
strongly depends on its radial profile, as already noticed by 
Nowak & Wagoner (1995). In particular, the result of this as- 
sumption is such that the models with stable and unstable disks 
in GBHs, as well as the model for the AGN disk, differ sub- 
stantially from each other even for the same prescriptions for 
the mass evaporation rate. This is because the disk thickness 
increases locally in the unstable zone of the radiation pressure 
dominated disk in GBH (cf. Fig[T3]). In AGN gas pressure dom- 
inated models the disk thickness in the outer rises slowly in the 
outer parts (cf. Fig[T4b. The slopes of the models (16) and (22) 
for a GBH and an AGN as given in Table|2]are both similar but 
this is the coincidence. If we calculate the AGN slope in the fre- 
quency range of log/ between -8.8 and -7.8, we obtain much 
flatter slope of 0.43. When we assume geometrically fixed dy- 
namo cells (b or c) the differences between the models remain. 
In particular, the conclusion based on analytical models that 
for Cb Pgas models the slope is 1 .0 both for GBH and AGN is 
not supported by the numerical computations: the slope is 0.76 
(model 26) and 0.93 (model 25), but both spectra show certain 
curvature. 

The prescription "B" for the mass evaporation rate gives al- 
ways much too strong emission in the corona (e.g. for 10 
black hole logLdisk ~ 38.5, and logLcor ~ 39), which in these 
models exceeds the disk luminosity by a factor of ~ 3. This 
is not consistent with the Soft State spectra being dominated 
by the disk component. Strong outflow from the disk to the 
corona in such models leads to the corona being optically thick 
for scattering, while the disk surface density drops consider- 
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Fig. 14. The ratio of the disk thickness to the radius as a func- 
tion of radius in AGN (Pgas heating). The black hole mass is 
10*^ Mo and the accretion rate is m = 0.35. 



ably in comparison with prescription "A" models for the same 
accretion rate, although it is never completely evaporated. We 
conclude that the assumption of the light speed instead of the 
Alfven velocity in energy flux transported to the corona does 
not lead to acceptable models. 

The evaporations laws "A" and "C" give disk-dominated 
models, and "C" (corresponding to modified inner boundary 
condition) gives always much flatter power spectra. Within the 
frame of these models, the magnetic dynamo parameter j6s 
plays also an important role. In the case of a stable disk in 
GBH and in AGN, it usually (models with f gas heating) has 
to be taken equal to 1.0, so that the dynamo can produce the 
substantial amplitude of X-ray (coronal) variability, supported 
by the observations of the Soft State sources (e.g. Treves et al. 
[TW Nowak et al. 'SOOT, Churazov et al. UMT\ Axelsson et 
al. 12006, for GBH; Markowitz |2003 Uttley & McHardylMll 
Shemmer et al.EESfor AGN). 

The /3s - I means that the global (large scale) magnetic 
field forms effectively and supports the energy transport be- 
tween the disk and corona. In the unstable GBH disks, the 
stochastic fluctuations of coronal emission are large even for 
the values of j6s smaller than 1.0, i.e. we can assume it to be 
0.5; otherwise, the corona is too strong. Note, that we are using 
here a constant value of the j6s parameter, but other possibili- 
ties have also been recently explored. For instance, Czerny et 
al. ( [2003] ) and Tout & Pringle ( (19921 ) suggested that Ps is pro- 
portional to H/r. Such a prescription steepens or flattens the 
PSD, depending on the radial profile of H/r. 

Finally, we stress that the mass of the central black hole also 
plays the role, and the PSD spectra corresponding to the ac- 
creting supermassive black holes are usually steeper in full nu- 
merical simulations of disk-corona evolution, in comparison to 



the analytical estimations. This is because the disk properties, 
when calculated properly, do depend on the black hole mass in 
a complex way. 

Analytically, the evaporation rate is proportional to BIvam 
(see Sect. 12.31 ). which, in the stable accretion disk, depends on 
the radius as r"^, regardless of the central mass. However, in 
the numerical simulation, when the coupling with the corona is 
included, this relation is modified. We find that for GBH disk 
the (averaged, i.e. not including the random changes) slope of 
mz(r) is ~ -3.28, while for AGN it is ~ -2.89 in the main 
inner part of the disk (up to log r ~ 15.8 cm for Mbh = IO^Mq; 
cf. Fig. [T4l i. Moreover, in the outer disk part (which is much 
more geometrically thick), the slope of this function changes 
dramatically and becomes positive: ~ 4.3. These properties are 
the reason for the systematically steeper PSD slopes in AGN, 
since in those disks the longer timescales are more pronounced 
than in Galactic sources. 

King et al. (2004) did not find any dependence on black 
hole mass, as these authors assumed that all the quantities scale 
with the BH mass in a self-similar way and the ratio H/r is 
constant. Mayer & Pringle (2006) did not study a difference be- 
tween the PSD spectra in stellar mass and supermassive black 
holes, and the dynamo and power spectra were calculated only 
for a fixed mass of 10 Mq. These authors discuss, however, 
the systematic differences in the instability picture. They find 
that for higher mass the disk becomes more unstable due to the 
larger range of accretion rates which are in the unstable regime. 

4.2. Observed PSD of the sources in the Soft States 

A few Galactic sources are permanently in the Soft State, many 
more sources (several transient X-ray sources, in particular) 
exhibit temporary Soft States (for a review, see Remillard & 
McClintock|2006|. 

The shape of the PSD was studied for a few sources in the 
Soft State (see Table[3ll. In most of the sources the power spec- 
trum is broad, with a slope 0.9 - 1.3. The power spectrum of 
XTE J 1650-500 is an exception since in this case the dominat- 
ing contribution to the power comes from a relatively narrow 
range of frequencies, around and below ~ 2 Hz. 

The most detailed analysis of the Soft State was done 
for Cyg X-1 by Axelsson et al. (120061 ). They analyzed 614 
power spectra by decomposition into two Lorentzians and a 
power law component with high frequency cut-off P{f) - 
af^^"''^ ^^V{-f I fc)- The distribution of the power law indices 
a for 523 successful fits is shown in Fig. [15] The mean slope 
IS, p - 0.92, with the dispersion of 0.18, and the mean high 
energy cut-off is 10.0 Hz, with the dispersion 9.4 Hz. 

Among AGN, it is widely believed that Narrow Line 
Seyfert 1 galaxies are high mass analogs of the GBHs in the 
Soft State (e.g. Pounds et al. [T995l McHardy et al. l2003l Uttley 
& McHardv [2005l Czernv l2003] ). In Table|4]we list objects with 
determined PSD slopes. In this case the power spectrum is not 
determined so accurately as in the case of Galactic sources, but 
the overall trend seems to be the same. In particular, the slope 
is about 0.8 - 1.3 both in the Soft State GBH and in NLSls, 
although better measurements is clearly needed. 
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Table 3. Power spectra of Galactic Black Holes in Soft State 



object 


low frequency slope 


high frequency slope 


frequency break 


reference 


Cyg X-1 


0.92 




~ 10 Hz 


Axelsson et al.|2006l 


XTE J 1650-500 


~ 0.06 


~ 1.4 


~ 2Hz 


Done & Gielinski 120051 


GRS 1915+105 


1.3 ±0.1 




> 10 Hz 


Rao etal.|2000| 


GRS 1915-1-105 


1.3 + 0.1 


1.61 + 0.02 


1 Hz 


Rao et al. 20001 


GRS 1915-1-105 


1.25 + 0.02 




> 10 Hz 


Rao etal. 120001 


LMC X-1 


~ 1 




0.2 Hz 


Nowak etal.120011 



Table 4. Power spectra of Narrow Line Seyfert 1 galaxies 



object low frequency slope high frequency slope frequency break reference 



MCG -6-30-15 


n s+u iEi 

"•°-0.4 


1 98+""" 


6.01^" X 10- 


^ Hz 


McHardy et al. 2005| 


NGC 4051 


1.1 (fixed) 




8.3+^' X 10- 


* Hz 


McHardy et al. 2004] 


Mrk 766 


a 1 


^•°-0.4 


a 5 X 10-* 


Hz 


Vaughan & Fabian 2003 1 


Ark 564 


1 94+0.03 


1 90+"" 


2.3t« X 10- 


3 Hz 


Papadakis etal., 20021 


NGC 5506 


~ 1.0 


> 1.7 


1.3!^;^ X 10- 


5 Hz 


Uttlev & McHardvl2005l 




0.6 0.8 1 1.2 1.4 



Fig. 15. The distribution of the slopes of the power law compo- 
nent in the power spectrum of Cyg X-1 in the Soft State (data 
taken from table 1 of Axelsson et al. 120061) 

Several models considered in our paper satisfy this re- 
quirement. Among the semi-analytical models, examples are: 
Ab-rad and Ab-gas, as well as Aa-gas and Ca-gas, which are 
roughly consistent with this requirement. Also, models Cb-rad 
and Cb-gas certainly provide good representation of the data. 
As for the full numerical models, several small mass models 
are acceptable (01, 03, 06, 07, 09, 13, 16 and 28 in Tabled, 
and four of the high mass models is satisfactory (19, 21, 22, 
25). 

5. Conclusions 

In this article we present the calculations of the time-dependent 
accretion disk plus corona model, incorporating the magnetic 



dynamo as a mechanism for the origin of the short-term vari- 
ability in hard X-rays. The current model is suitable for the 
Soft State X-ray binaries, as well as their supermassive black 
hole analogues: Narrow Line Seyfert 1 galaxies. The radiation 
spectrum is dominated by the thermal soft component, i.e. the 
accretion disk, which may either be constant in time, or vari- 
able on relatively long timescale due to the instability induced 
by the radiation pressure domination. The hard spectral com- 
ponent caiTies less energy, and is produced by the accreting 
corona above the disk. The coronal emission is variable on 
short timescale due to the the rapidly changing magnetic field 
in the accretion disk. The key element of the model is there- 
fore the coupling between the disk and corona, which due to 
the mass evaporation driven by the magnetic flux allows the 
perturbations generated locally in the disk to be transported to 
the corona. Since the viscous timescale in the corona is much 
shorter than in the disk, these perturbations are not smeared and 
a broad power spectrum is produced. 

By studying our model both analytically and numerically, 
as well as by the comparison with the available observational 
data, we drive the following conclusions: 

1 . The radiation pressure instability in the underlying accre- 
tion disk leads to steeper power density spectrum of its fluc- 
tuations, as well as to the larger amplitudes. 

2. The assumption of the light speed instead of the Alfven ve- 
locity in energy flux transported to the corona does not lead 
to acceptable models. 

3. The PSD spectra never look like a broken power law, they 
turn off slowly at higher frequencies and frequently show 
some curvature at low frequencies as well. 

4. The PSD spectra corresponding to the accreting supermas- 
sive black holes for most models are steeper than for the 
galactic black holes. This result has been obtained in full 
numerical simulations of disk-corona evolution, while the 
analytical estimations give PSD independent on BH mass. 

The latter point might seem inconsistent with the AGN data 
obtained so far, since their PSD spectra seem to have slopes 
similar to these of Galactic sources. However, one should take 
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into account the still poor statistics and large errors of these ob- 
servations. We conclude that further verification of the models 
requires more detailed observations and analysis of the long- 
term variability in AGN. 

Acknowledgements. We thank to Piotr Zycki for helpful discussions. 
Part of this work was supported by grant 1P03D00829 of the Polish 
State Committee for Scientific Research. 

References 

Abramowicz, M.A., Lanza, A., Spiegel, E.A., Szuszkiewicz, E., 1992, 
Nat., 356, 41 

Abramowicz, M.A., Czerny, B., Lasota, J. -P., Szuszkiewicz, E., 1988, 

ApJ, 332, 646 
Abrassart, A., Czerny, B., 2000, A&A, 356, 475 
Axelsson, M., Borgonovo, L., Larsson, S., 2006, A&A, 452, 975 
Balbus, S.A., Hawley, J.R, 1998, RvMP 70, 1 
Bath,, G.T., Evans, W.D., Papaloizou, J., 1974, MNRAS, 167, 7P 
Bisnovatyi-Kogan, G.S., Blinnikov, S.I., 1977, A&A, 59, 111 
Boiler Th., Brandt W.N., Fabian A.C., Fink H.H., 1997, MNRAS 289, 

393 

Canizares, C.R., Oda, M., 1977, ApJL, 214, LI 19 

Churazov, E., Gilfanov, M., Revnivtsev, M., 2001, MNRAS, 321, 759 

Czerny, B., 2003, in Active Galactic Nuclei: from Central Engine to 

Host Galaxy, Eds.: S. Collin, F. Combes and I. Shlosman. ASP 

Conference Series, Vol. 290, p. 59. 
Czerny, B., Nikolajuk, M., Rozanska, A., Dumont, A.-M., Loska, Z., 

Zycki. PT, 2003, A&A, 412, 317 
Deufel, B., Spruit, H.C., 2000, A&A, 362, 1 
Done, C, Gierlinski, M., 2005, MNRAS, 364, 208 
Galeev, A.A., Rosner, R., Vaiana, G.S., 1979, ApJ, 229, 318 
Gleissner, T, Wilms, J., Pottschmidt, K., Uttley, P., Nowak, M.A., 

Staubert, R., 2004, A&A, 414, 1091 
Goosmann, R.W. et al., 2006, astro-ph/0604T56l 
Hawley, J. F, Gammie, C. F, Balbus, S. A., 1995, ApJ, 440, 742 
Janiuk, A., Czerny, B., 2005, MNRAS, 356, 205 
Janiuk, A., Czerny, B., Siemiginowska, A., 2000, MNRAS, 318, 180 
Janiuk, A., Czerny, B., Siemiginowska, A., 2002, ApJ, 576, 908 
Janiuk, A., Zycki, PT, Czerny, B., 2000, MNRAS, 314, 364 
Karas, V., Czerny, B., Abrassart, A., Abramowicz, M.A., 2000, 

MNRAS, 318, 547 
King, A.R., Pringle, J.E., West, R.G., Livio, M., 2004, MNRAS, 348, 

111 

Komissarov, S.S., Barkov, M., Lyutikov, M., 2006,|astro-ph/0606375| 
Kotov, O., Churazov, E., Gilfanov, M., 2001, MNRAS, 3277 799 
Kuperus, M., lonson, J.A., 1985, A&A, 148, 309 
Lawrence, A., Watson, M.G., Pounds, K.A., Elvis, M., 1987, Nat., 
325, 694 

Livio, M., in Wickramasinghe D. T., Bicknell G. V., Ferrario L., 
eds, ASP Conf. Ser. Vol. 121, Accretion Phenomena and Related 
Outflows, lAU Colloq. 163. Astron. Soc. Pac, San Francisco, p. 
845. 

Lovelace, R.V.E., Romanova, M.M., Newman, W.I., 1994, ApJ, 437, 
136 

Lyubarskii, Yu. E. 1997, MNRAS, 292, 679 

Markowitz, A., et al., 2003, ApJ, 593, 96 

Mayer, M., Pringle, J.E., 2006, MNRAS, 368, 379 

McHardy, I.M., Czerny, B., 1987, Nat., 325, 696 

McHardy, 1. M., Papadakis, I. E., Uttley, P., Page, M. J., Mason, K. O. 

2004, MNRAS, 348, 783 
McHardy, I.M., Gunn, K.F, Uttley, P, Goad, M.R., 2005, MNRAS, 

359, 1469 



McHardy, I.M., Mason, K., Page, M., 2003, in van den Heuvel E. P. 

J., in 't Zand J. J. M., Wijers R. A. M. J., eds. The Restless High 

Energy Universe, Elsevier 
McKeman, B., Yaqoob, T, 1998, ApJ 501, L29 
Merloni, A., 2003, MNRAS, 341, 1051 
Meyer, F, Meyer-Hofmesiter, E., 1994, A&A, 288, 175 
Meyer-Hofmesiter, E., Meyer, F, 2001, A&A, 380, 739 
Meyer-Hofmeister, E., Liu, B.-F, Meyer, F, 2005, A&A, 432, 181 
Misra, R., Taam, R.E., 2002, ApJ, 573, 764 
Muchotrzeb, B., Paczynski, B., 1982, AcA, 32, 1 
Mushotzky, R.F, Done, C, Pounds, K.A., 1993, ARA&A, 31, 717 
Nayakshin, S., Rappaport, S., Melia, F, 2000, ApJ, 535, 798 
Negoro, H., Kitamoto, S., Takeuchi, M., Mineshige, S., 1995, ApJ, 

452, L49 

Nowak, M.A., Wilms, J., Heindl, W.A., Pottschmidt, K., Dove, J.B., 

Begelman, M.C., 2001, MNRAS, 320, 316 
Nowak, M. A., Wagoner, R. V., 1995, MNRAS, 274, 37 
Paczynski, B., Bisnovatyi-Kogan, G., 1981, AcA, 31, 283 
Paczynski, B., Wiita, PL, 1980, A&A, 88, 23 

Papadakis, I.E., Brinkmann, W., Negoro, H., Gliozzi, M., 2002, A&A, 
382, LI 

Pottschmidt, K., Wilms, J., Nowak, M. A., et al. 2003, A&A, 407, 
1039 

Pottschmidt, K., Chernyakova, M., Zdzziarski, A.A., Lubinski, P., 

Smith, D.M., Bezayifl^, N., 2006, A&A, 452, 285 
Pounds, K.A., Done, C, Osborne, J.P, 1995, MNRAS, 277, L5 
Poutanen, J., Fabian, A., 1999, MNRAS, 306, L31 
Pringle, J.E., Rees, M.J., Pacholczyk, A.G., 1973, A&A, 29, 179 
Rao, A.R., Yadav, J.S., Paul, B., 2000, ApJ, 544, 443 
Remillard, R.A., McClintock, J.E., 2006, astro-ph/0606352 
Revnivtsev, M., Gilfanov, M., Churazov, E., 2000, A&A, 363, 1013 
Rozanska, A., Czerny, B., 2000, A&A, 316, 1055 
Shakura, N.I., Sunyaev, R.A., 1973, A&A, 24, 337 
Shakura, N.L, Sunyaev, R.A., 1976, MNRAS, 175, 613 
Shemmer, O., Brandt, W.N., Vignali, C, Schneider, D.P, Fan, X., 

Richards, G.T., Strauss, M.A., 2005, ApJ, 630, 729 
Smith, D. M.; Heindl, W. A.; Swank, J. H., 2002, ApJ, 569, 362 
Stone, J.M., Hawley, J.F, Gammie, C.F, Balbus, S.A., 1996, ApJ, 463, 
656 

Szuszkiewicz, E., Miller, J.C., 1997, MNRAS, 287, 165 

Szuszkiewicz, E., Miller, J.C., 1998, MNRAS, 298, 888 

Taam R.E., Lin D.N.C., 1984, ApJ, 287, 761 

Terrel, N.J., Jr, 1972, ApJL, 174, L35 

Tout, C.A., Pringle, J.E., 1992, MNRAS, 259, 604 

Tout, C. A., Pringle, J. E., 1996, MNRAS, 281, 219 

Treves, A., Belloni, T., Chiapetti, L., Maraschi, L., Stella, L., Tanzi, 

E.G., van der Khs, M., 1988, ApJ, 325, 119 
Uttley, P, McHardy, I.M., 2001, MNRAS, 323, 26 
Uttley, P, McHardy, I.M., 2005, MNRAS, 363, 586 
van der Klis, M., 1989, ARA&A, 27, 517 
Vaughan, S., Fabian, A.C., 2003, MNRAS, 341, 496 
Witt, H.J., Czerny, B., Zycki, PT, 1997, MNRAS, 286, 848 
Zycki, PT, CoUin-Souffrin, S., Czerny, B., 1995, MNRAS, 277, 70 

Appendix 

In the computations of the time evolution of the disk corona 
system, the energy equation (Eq. 23) is simphfied due to the 
vertical integration over the disk height (see e.g. Taam & Lin 
1984). This approach neglects the components connected with 
the time evolution of the disk thickness, which means an as- 
sumption that the density decreases significantly at the disk 
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surface. Here we derive a full energy equation that includes 
other the necessary terms, e.g. assumes that the vertical distri- 
bution of density in the disk is uniform. Then we check if our 
simplified treatment is justified. 

The energy equation in its basic for is: 



dS 

T — S = - G- 



(25) 



where S is the entropy and other terms have the same meanings 
as above. 

From the thermodynamical relation we have: 



TdS ^-[{12- \0.5l3)dlnT - (4 - 3/3)<i/np] (26) 
P 

and therefore 

PS dliiT dlnp 

—[(12 - 10.5;6)— - - (4 - 3^)— f ] = - Q . (27) 
p dt dt 



Dividing the above equation by P2/p/(12 - lQ.5/3) we have: 

(28) 



dlnT _ 4-3/3 dlnp ^ (g^ - Q-)p 



dt 12-10.5^ dt PZ(12- 10.5/?)' 
and substituting p with 'L/H, this turns into: 
dlnT _ 4-3/3 dlnp _^ (g^ - Q ) 



dt 12-10.5/3 dt PHiU - 10.5/3) 



(29) 



Since we use dlnT - dT/T, the above equation is divided 
by T. However, using dlnp = t/lnS is an simplification that 
we have done in the numerical calculations presented above. 
To be strict, we should take rather: 



P = 



H 



which is equivalent to: 
Inp = In 2 - In // 
and 

dlnp = dlnY, - dlnH 



(30) 



(31) 



(32) 



In order to have dlnH as a function of T and S, we take the 
hydrostatic balance equation: 



P 



and 

P = Q^//V = 



H 



(33) 



(34) 



The pressure is a sum of gas and radiation pressures, and 
therefore: 



k 1 , k I.T 1 , 

P = pT + -aT^ = + -aT^ 

nip/u 3 nipij H 3 

From the two equations above, we have: 

k XT 1 , , 



Mpid H 3 



(35) 



(36) 



Now we differentiate our equation, taking into account that 
the first term above is equal to /3P, and the second term is ( 1 - 
/3)P. Starting with: 

k ^dl. 
nipju H 



-T"-^ -r^XT% + z^X^ + UT^dT 



m^fi H 3 



d(Q})HI. + Q}dHi: + a^Hdl. 
and after a couple of arithmetical operations we have: 



(37) 



(/?- l)£/lnS-H(4-3^)t/lnr-£/ln(Q^) = {l+/3)dlnH (38) 

From the above equation we obtain dlnH and put it into the Eq. 
[32l which gives: 



dlnp 



4-3/3 



1 



dlnl.-- -dlnT + - dln{n^). (39) 



1+/3 1+/3 1+/3 

This relation is now substituted into the equation |27] 

PH[(12- 10.5/3)dlnT -^[2dlnX - (4- - 3p)dlnT + 

flf/n(Q2)]] = e+ - e_ (40) 
and after a few transformations we have: 

dlnT = '-^^dlnX.'-^dln^^).^l^^m:^,,^ 
Ab Ab PHAb 

where: 

Ab = [(12 - 10.5yS)(l + /3) + (4 - 3,6)2] ^^2) 

Finally, we delogarythimize the equation: 
dT = 2r(4 - 3/3)^^^7X4-3^ r(l . /3XQ^ - Q,) 

AbS Ab PHAb 

and substitute the full time derivatives of dT and dl, with partial 
ones: 

^r = f.v,f: (44) 
ot ar 

dl. dl 

fl^S = — + V,— (45) 
at Or 

This equation is more complicated than the one used in our 
numerical calculations. Therefore, let us now check, how the 
numerical results depend on the adopted formula for dT/dt. 
We performed several test runs with the above relation and we 
found that the results only slightly differ from the ones with 
the simplified formula used throughout this work. The test runs 
were performed for models 01, 04, 07 and 19 as in the Table |2l 
In general, for the models in which the disk was stable (07 and 
19), the value of its constant luminosity was exactly the same 
in the test runs as in the Table. For the models in which the disk 
was unstable due to the radiation pressure (01 and 04), the am- 
plitudes of the luminosity outbursts and their durations changed 
by about 0. 1 %. The same range of change occurred for the am- 
plitudes of the coronal variability in all of the tested models. 
For the slopes of the power density spectra, they changed by 
9%, 3%, 2% and 8%, for the test models 01, 04, 07 and 19, 
respectively. These differences were most probably due to the 
numerical reasons (the time step used in our model is variable, 
depending on the rate of change of temperature and density in 
the disk). We conclude, that the approximation used in case of 
the energy equation in the time dependent numerical calcula- 
tions is viable and does not influence our results in a great deal. 



